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Abstract 


y(t)  = 6(t')  = x1(t')  . (4) 


This  paper  deals  with  the  mathematical  modeling  of  a 
dc  motor  and  quadratic  optimal  regulator  system. 
MATLAB  based  program  is  used  to  design  the 
quadratic  optimal  regulator  system  for  the  dc  motor. 
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mathematical  modeling,  MATLAB. 

1.  Introduction 

We  know  that  it  is  possible  to  design  state  feedback 
controllers  using  using  eigenvalue  (pole)  placement 
algorithms.  For  single  input  systems,  given  a set  of 
desired  eigenvalues,  the  feedback  gain  to  achieve  this  is 
unique  (as  long  as  the  system  is  controllable).  How  does 
one  utilize  this  freedom?  A more  fundamental  issue  is 
that  the  choice  of  eigenvalues  is  not  obvious.  For 
example,  there  are  tradeoffs  between  robustness, 
performance,  and  control  effort. 

Linear  quadratic  (LQ)  optimal  control  can  be  used  to 
resolve  some  of  these  issues,  by  not  specifying  exactly 
where  the  closed  loop  eigenvalues  should  be  directly,  but 
instead  by  specifying  some  kind  of  performance 
objective  function  to  be  optimized. 

The  system  that  minimizes  the  selected  performance 
index  is,  by  definition  optimal.  Although  the  controller 
may  have  nothing  to  do  with  “optimality”  in  many 
practical  applications,  the  important  point  is  that  the 
design  based  on  the  quadratic  performance  index  yields  a 
stable  control  system.[5],[6],[7]. 


We  chose  the  different  state  variable  of  the  DC  Motor 
asx1(t)  = 0(t),  x2(t)  = a )(t)  = de^  and  x3(t)  = 

ia(t). 

Where  ea(t)  = armature  voltage,  ia  (t)  = armature 
current,  0(t)  = motor  shaft  angle,  ~~  = co(t)=  shaft 
speed,  J = moment  of  inertia  of  the  rotor,  B = viscous 
frictional  constant,  Lm=  inductance  of  armature 
windings,  Rm=  armature  winding  resistance,  Kt=  motor 
torque  constant,  Km=  motor  constant. 

Here  the  motor  speed  co  (t)  is  controlled  by  varying  the 
armature  voltage  ea(t).  Hence  ea(t)  is  the  input  variable 
and  co  (f)  is  the  output  variable. 


Hence  the  state  model  of  DC  Motor  is  derived  from 
equations  (1),  (2),  (3)  and  (4)  as  follows 
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Let,  the  motor  parameters  (coefficient  of  differential 
equations)  are  assigned  to  be  Lm  = 0.5  H,  Kt  = 0.01  N- 
m/A,  Km=  0.01  V-sec/rad,  J = 0.01  kg-m2,  B = 0.1  N-m- 
sec/rad,  Rm=  1 Q. 


This  paper  is  mainly  discussing  about  the  modeling  of  dc 
motor  [1],  [2],  and  Quadratic  Optimal  Regulator. 

MATLAB  based  analysis  is  also  shown  at  the  end. 

2.  Mathematical  Modeling  of  DC  Motor 


The  different  equations  related  to  DC  Motor  are  given 
below: 
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Thus  the  state  model  of  dc  motor  is  derived  using  motor 
parameters  and  equation  (5)  and  (6)  as  follows: 
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3.  Mathematical  modeling  of  Quadratic 
Optimal  Regulator 

The  given  system  equation  is 

x = Ax  + Bu  (9) 

This  determines  the  matrix  K of  the  optimal  control 
vector 

u = -Kx  (10) 

So  as  to  minimize  the  performance  index 
J = /“(JTQX  + uT Ru)  dt  (11) 

Where  Q is  a positive  definite  (positive  semi  definite) 
Hermitian  or  real  symmetric  matrix  and  R is  a positive - 
definite  Hermitian  or  real  symmetric  matrix. 

It  can  be  proved  that  if  A-BK  is  a stable  matrix,  there 
exists  a posite  definite  matrix  P that  satisfies  the 
following  equation: 

(A  - BK)TP+P(A-BK)=-(Q+KTRK)  (12) 

Since  all  eigen  values  of  A-BK  are  assumed  to  have 

negative  real  parts,  we  have  x(oo)  ->  0,  therefore  we  can 
write 

J = XT(0)PX(0)  (13) 

Thus  performance  index  J can  be  obtained  in  terms  of  the 
initial  condition  X (0)  and  P. 

Thus  the  optimal  matrix  K is  given  by 
K = R~1BtP  (14) 

The  matrix  in  equation  (14)  must  satisfy  equation  (12)  or 
the  following  equation: 

AtP+PA-PBR~1BtP+Q=Q  (15) 


K = 


1.0000 

0.1069 

0.4515 

14.6323 

1.4515 

0.5000 

1.4515 

0.1945 

0.0535 

0.5000 

0.0535 

0.2258 

E = 

-0.0705 

-2.8373 

-9.9952 

5.  Conclusion 

The  result  has  shown  the  gain  matrix  K,  eigenvalue 
vector  E,  and  matrix  P of  the  quadratic  optimal  regulator 
system  for  dc  motor  because  matrix  A-BK  is  a stable 
matrix.  It  is  important  to  note  that  for  certain  systems 
matrix  A-BK  cannot  be  made  a stable  matrix,  whatever  K 
is  chosen.  In  such  a case,  there  does  not  exist  a positive- 
definite  matrix  P for  the  matrix  Riccati  equation. 


The  above  equation  is  called  the  reduced  -matrix  Riccati 
equation. 

The  design  steps  may  be  stated  as  follows: 

a)  Solve  equation  (15),  the  reduced  Riccati  equation  for 
matrix  P. 

b)  Substititute  this  matrix  P into  equation  (14).  The 
resulting  matrix  K is  the  optimal  matrix. 

4.  MATLAB  Based  Quadratic  Optimal 
Regulator  for  dc  motor 

The  given  MATLAB  command  is  used  to  solve  the 
continuous  time,  linear,  quadratic  regulator  problem  and 
the  associated  Riccati  equation. 

» A=  [0  10:0 -10  1;0  -0.02  -2]; 

B = [0;0;2]; 

Q=  [1  0 0:0  1 0;0  0 1]; 

R = [ 1]; 
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